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We present a search for the possible I(J P ) = 0(2 + ) tetraquark state with ssss quark content in 
quenched improved anisotropic lattice QCD. Using various local and non-local interpolating fields we 
determine the energies of ground-state and second ground state using variational method. The state 
is found to be consistent with two-particle scattering state, which is checked to exhibit the expected 
volume dependence of the spectral weights. In the physical limit, we obtain for the ground state, 
a mass of 2123(33) (58) MeV which is higher than the mass of experimentally observed /(2010). 
The lattice resonance signal obtained in the physical region does not support a localized J p = 2 + 
tetraquark state in the pion mass region of 300 — 800 MeV. We conclude that the 4q system in 
question appears as a two-particle scattering state in the quark mass region explored here. 

PACS numbers: 11.15.Ha, 11.30.Rd,12.38.Ge 
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I. INTRODUCTION 

The concept of multi-quark hadrons has received re- 
vived interest due to the narrow resonances in the spec- 
trum of states. Recently, several new particles were ex- 
perimentally discovered and confirmed as the candidates 
of multi-quark states. These discoveries are expected 
to reveal new aspects of hadron physics. Among these 
discoveries, the tetra-quark systems are also interesting 
in terms of their rich phenomenology, in particular for 
mesons which still remain a most fascinating subject 
of research. The 4q states are interesting in terms of 
the recent experimental discoveries of X(3872) [l], 0,3 , 
Y(4260) [| and L> s (2317) @, @], which are expected to 
be tetra-quark candidates. 

The Particle Data Group lists 2 tensor mesons with 
masses in the range 1.9 — 2.2 GeV/c 2 and considers them 
as well-established. The 2++ candidates, / 2 (1950) @,S, 
HI and /2(2010) [l(| are isosinglet. The relevant channels 
of decay are KK and r\r\ for the /2(1950) and <j><j> and KK 
for /2(2010). Due to their KK decay, one would expect 
/2(1950) and /2(2010) are very likely one state; the mass 
shift could be a measurement error or could be caused 
by the KK threshold. However, the results for / 2 (2010) 
favour an intrinsically narrower state, strongly coupled to 
4>4> and weakly coupled to the other channel for allowed 
s-wave decays. Following the recent re-analysis of the 
BNL data [n[ we discuss the state /2(2010) as a s 2 s 2 
state. 

The multi-quark states have been investigated in lat- 
tice QCD studied with somewhat mixed results [HI, 0, 
Hil . [ToT ]. At the present status of approximations, lat- 
tice QCD seems to provide a trustworthy guide into un- 
known territor y in tetra-quark hadron physics [l6l . UtI 
EE EI HE HUM HE [H, HE HH. Using the quenched 
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approximation, and discarding quark-antiquark annihila- 
tion diagrams, we construct s 2 s 2 sources from multiple 
operators. Note that we are working in quenched approx- 
imation which in principal is unphysical. However, previ- 
ous lattice results on masses and decay constants turn out 
to be in good agreement with experimental values [27| . 
This seems to suggest that it is plausible to use quenched 
lattice QCD to investigate the mass spectra. We exclude 
the processes that mix qq and q 2 q 2 and allow the quark 
masses to vary from small to large values. In the absence 
of quark annihilation, we do not expect any mixing of 
q 2 q 2 with pure glue. Thus we can express the q 2 q 2 corre- 
lation functions in terms of a basis determined by quark 
exchange diagrams only (ignoring the single, double and 
annihilation diagrams among Wick's contractions). An- 
other important question is whether the interpolating op- 
erator one uses has a significant overlap with the state in 
question. To construct an interpolating field which has 
significant overlap with the 4g system, we adopt the so- 
called variational method to compute 2x2 correlation 
matrix from two different interpolating fields and from 
its eigenvalues we extract the masses. Thus, assuming 
that the quenching uncertainties do not effect our con- 
clusions dramatically, we investigate the optimized corre- 
lation function and use it to examine lowest-lying tetra- 
quark resonance as f2(ssss) states in the spectrum of 
2x2 correlation matrix. 



II. LATTICE STUDY FOR THE f 2 (ssss) 

The simplest local interpolators can be written in 
terms of colour-singlet configuration of a product of 
colour-neutral meson interpolation fields. We propose 
a non.-(fxj) interpolating field to extract the fi(ssss) 
tetraquark state. This choice is designed to maximize 
the possibility to observe attraction between tetraquark 
constituents at relatively large quark masses. With a 
4><p operator it is possible that there is a small amount of 
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the compact Aq component in the two-body interpolating 
field since the interpolator may contain a large contam- 
ination of 4><P scattering states. We adopt the simplest 
non-</x/)- type interpolator of the form 



Oi(x) = - 



q <-> Q,q <-> Q 



(1) 



with spin I(J P ) = 0(2+) . In the nonrelativistic limit the 
above non-two state particle can not be decomposed into 
</>(/). Thus the 4g state can be singled out as much as pos- 
sible and the results are less biased by the contamination 
of two-state scattering states. 

The other type of interpolating field is one in which 
quarks and anti-quarks are coupled into a set of diquark 
and antidiquark, respectively and has the form 



2 {x) = e abc [q^CTQc] e ade [q d CTQ^ 



(2) 



Accounting for both colour and flavour antisymmetry, 
possible Ts are restricted within 75 and For T = 
(i = 1) 2, 3), the above diquark operator transforms 
like J p = 1~. For concreteness, we simulate the flavour 
combination [ss] and [ss]. 

To extract energies E n of s 2 s 2 We compute the 2x2 
correlation matrix 



C i ,(t) = (^tr[((O i )(x,t)O i )(0,0)) / 



(3) 



where the trace sums over the Dirac space, and the sub- 
scripts / and U denote fermionic average and gauge fi eld 
ensemble average, respectively. Following [T2l.fl6l.l28j we 
solve the eigenvalue equation 



C(toH(*o) = Xk(t )v k {t ) 



(4) 



to determine the eigenvectors Vk(to). We use these eigen- 
vectors to project the correlation matrices to the space 
corresponding to the n largest eigenvalues X n {to) 



Cm) = (v i ,C(t)v j ), 



h3 



(5) 



and solve the generalised eigenvalue problem equation 
for the projected correlation matrix Cfj. The resulting 
large-time dependence of the eigenvalues X n (t) allows a 
determination of ground and excited-state energies. The 
mass can be extracted by a hyperbolic-cosine fit to A n (t) 
for the range of t in which effective mass 



M eff (t)= In 



A(t) 



A(i + 1) 



(6) 



attains a plateau. In order to show the existence or ab- 
sence of the signature of tetraquark resonance on lattice, 
we establish lowest and as well as the second-lowest en- 
ergy levels for our 4q system. 

Using a tadpole-improved anisotropic gluon action 
[29I ] . we generate quenched configurations on two lattice 



volumes 16 3 x 64 and 16 3 x 80 (with periodic boundary 
conditions in all directions). After discarding the initial 
sweeps, a total of 200 configurations are accumulated for 
measurements at (3 = 4.0 Quark propagators are com- 
puted by using a tadpole- improved clover quark action on 
the anisotropic lattice [3(| ■ All the coefficients in the ac- 
tion are evaluated from tree-level tadpole improvement. 

The bare mass of the strange quark is determined 
by extracting the mass of the vector meson M^. At 
m q at = 0.066, we obtain n t = 0.2404, which produces 
a mass for the <f> of 1.237(2) in lattice units. Using the 
mass of the nucleon in the chiral limit, we find that the 
ratio M<j,/M N at the chiral limit is 1.059 ± 0.014, which 
is in good agreement with the physical ratio of 1.087. 
This verifies that the strange bare quark mass of 0.07 
used is very close (within 3%) to the physical strange 
quark mass. The quark propagators are then computed 
at seven values of the hopping parameter k t which cover 
the strange quark mass region of m s < m q < 2m s , i.e., 
a t m q = 0.07,0.075,0.08,0.09,0.105,0.115,0.12. Inspired 
by the good agreement of the ratio with the experimental 
value, the scale was set alternatively by M^/Mjy- Using 
the experimental value 938 MeV for the nucleon mass, 
the spacing of our lattice is a s = 0.473(2) fm. 



III. RESULTS AND DISCUSSION 

Figs. Q]illustrate the two lowest energy levels extracted 
by fitting the effective masses over appropriate t intervals. 
The ground state eigenvalues show a conventional time- 
dependence near t ~ T/2 and hence the mass can be 
accurately extracted using Eq. © . We choose one "best 
fit" which is insensitive to the fit range, has high con- 
fidence level and reasonable statistical errors. We then 
confirm this by looking at the plateau region of the cor- 
relator. Statistical errors of masses are estimated by the 
jackknifc method and the goodness of the fit is gauged by 
the x 2 /Ndf, chosen according to criteria that x 2 /Ndf 
is preferably close to 1.0. 

The effective mass is found to be stable using different 
values of t in Eq. which suggests that the ground 
state in question is correctly projected. Suppressing any 
data point which has error larger than its mean value, 
the possible plateau is seen in the region 5 < t < 12 
with reasonable errors, where the single-state dominance 
is expected to be achieved. Fitting the effective mass in 
the window t = 6 — 11 is found to optimize the x 2 /Ndf- 
To avoid the clutter in Fig. [T]we do not show the points 
at larger t values which have larger error bars, and have 
little or no influence on the fits. The best fit curve to 
the 4g data has x 2 / n df = 0.87. The results for the 
masses corresponding to the various values of the hopping 
parameter n t are tabulated in Table U 

To interpret the ground state in terms of signatures 
of a lattice resonance, we look at three possible sce- 
narios. First, we extract the mass splitting between 
the tetraquark 0(2 + ) and the noninteracting 4> + (j> two- 
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FIG. 1: Effective mass of the I (J) = 0(2+) colour -singlet 
lowest-lying ground state. The data correspond to ~ 361 
MeV (triangles) and 824 MeV (circles) . 



TABLE I: The masses of the 4q, Kaon and 
lattice units, for various values of k±. 



states, in the 





M 4q 


M 4 * 9 


M K 




0.2410 


2.187(4) 


2.223(6) 


0.579(2) 


1.029(3) 


0.2420 


2.055(7) 


2.083(9) 


0.507(5) 


0.973(6) 


0.2435 


1.951(13) 


1.973(11) 


0.475(8) 


0.926(11) 


0.2440 


1.812(19) 


1.825(21) 


0.440(5) 


0.860(15) 


0.2450 


1.724(24) 


1.751(37) 


0.417(13) 


0.822(28) 


0.2455 


1.655(27) 


1.673(44) 


0.401(14) 


0.787(17) 


0.2462 


1.591(23) 


1.606(49) 


0.383(19) 


0.759(24) 



particle state and compare our results to that derived in 
quenched chiral perturbation theory 1 . Fig. [5] shows the 



mass difference AM = M^ q — 2M^, together with the 
quenched one- loop energy shift in the finite box (3lj , as a 
function of m^L for the lowest 4g state from 16 3 x 64 lat- 
tice in our calculation. We obtain the results for one-loop 
energy shift by interpolating the coefficients A^m^L) 
and Boim-^L) listed in Ref. 31( for the range of m^L 
appropriate for our calculation on 16 3 x 64 lattice for 
8 = 0.12 and 0.15. 

We see clearly that the masses derived for the 
tctraquark state are consistently higher than the lowest 
two-particle state. The mass difference is over 100 MeV 
at small quark masses and weakly dependent on m^L. 
The positive mass difference observed in this range of 
pion mass suggests that the observed signal is unlikely 
to be a tetraquark. We also notice that our data are 
reasonably consistent with one-loop quenched perturba- 
tion results [3ll for m^L > 4.3 for 8 = 0.12 and 0.15. 




1 Since we are using the quenched approximation, the extraction 
of energy shift in a finite box using full QCD one-loop chiral 
perturbation theory is not applicable 



FIG. 2: The energy shift of the lowest I(J P ) = 0(2+) state 
as a function of m n L. The solid and dashed lines correspond 
quenched one-loop chiral perturbation results for 8 = 0.12 
and 0.15, respectively. 



It is interesting to note that our results are consistent 
with quenched one-loop results despite the fact that the 
disconnected contributions were not inserted in our cal- 
culation. This implies that disconnected correlator has 
very small or negligible contribution than the connected 
one at several time separations. 

To confirm or discard the signature observed in Fig. [21 
we examine the second scenario, i.e., the volume depen- 
dence of the spectral weight of these states. Theoreti- 
cally, if the state is a genuine resonance, then its spectral 
weight should be almost constant for any lattices with the 
same lattice spacing. On the other hand, if it is a two- 
particle scattering state, then its spectral weight has an 
explicit 1/V 3 dependence [IH. In the following, we shall 
use the ratio of the spectral weights on two spatial vol- 
umes 16 3 and 20 3 to discriminate whether the hadronic 
state in question is a resonance or a scattering state. 

Fig. [3] shows the ratio (R = W16/W20) of spectral 
weights of the lowest state and second-lowest state, ex- 
tracted from the time-correlation function of variational 
matrix as a function of m^. Since our two lattice sizes are 
16 3 x 64 and 20 3 x 80, the spectral weight ratio for a two- 
particle state should be Wi 6 /W 2 o = V20/V1Q = 1.95. Wc 
see that the ratio R for the lowest state clusters around 
1.0 for <G [0.5,0.8], which implies that there exists a 
2 + resonance with quark contents (ssss). 

On the other hand, for smaller quark masses, R begins 
to deviate from 1.0 with larger errors, suggesting that 
this state is a scattering state. Since none of our oper- 
ators has scalar meson component, the possibility that 
this might be due to quenching effects at smaller quark 
masses is highly unlikely. Thus one can safely ignore the 
possibility of R being consistent with 1.0 if one incor- 
porates internal quark loops with larger volumes. This 
type of flip-flop between the 4q state and the two-c/> state 
might be a flux -tube recombination between two 4> at 
some diquark and internal quark separations. This can 
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FIG. 3: Spectral weight ratio W16/W2 as a function of m% 
for the lowest state (solid circles) and next lowest state (solid 
triangles) . 



be verified by analyzing the Aq potential of the tetraquark 
system. We do not intend to pursue such an analysis here 
since this is not the focus of our present study. The spec- 
tral weight ratio of the first excited state turns out to be 
consistent with 1.95, confirming our speculation that it 
is two-particle scattering state. The two states are rea- 
sonably well separated compared to the decay width of 
/(2010). 

Finally, the mass differences extracted can be extrapo- 
lated to the physical limit, which is the next important is- 
sue [33j . Since quenched spectroscopy is quite reliable for 
mass ratio of stable particles, it is physically even more 
motivated to extrapolate mass ratios rather than masses 
or mass differences. This allows for the cancellation of 
systematic errors since the hadron states are generated 
from the same gauge field configurations and hence sys- 
tematic errors are strongly correlated. We use a set of 
data points with smallest m 2 to capture the chiral log be- 
haviour. Fig. U collects and displays the resulting mass 
ratios, illustrated in Tableim extrapolated to the physical 
limit using linear and quadratic fits in m 2 . . The difference 
between these two extrapolations gives some information 
about systematic uncertainties in the extrapolated quan- 
tities. Performing such extrapolations to mass ratios, we 
adopt the choice which shows the smoothest scaling ba- 
haviour for the final value, and use others to estimate the 
systematic errors. 

The data at smallest five quark masses behave almost 
linearly in m 2 and both the linear and quadratic fits es- 
sentially gave the identical results. The contributions 
from the uncertainties due to chiral logarithms in the 
physical limit are seen to be significantly less dominant. 
The mass difference AM is ~ 100 MeV at the smaller 
quark masses, and weakly dependent on m„. The sig- 
nature of repulsion at quark masses near the physical 
regime would imply no evidence of the resonance in the 
J = 2 channel. If this mass difference from two-</> thresh- 
old can be explained by the two-0 interaction, then the 
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FIG. 4: Extrapolation of the mass ratio AM/Mk for lowest 
state (open circle) and second lowest state (solid triangles) 
to the physical limit at a a = 0.473 fm. Also are shown the 
mass ratio Mk/M^ (solid diamonds). The dashed lines are 
the linear fits in m£ to the data. 



s 2 s 2 state can be regarded as a two-0 scattering state. 



TABLE II: Hadron mass ratios at various pion masses at 
a s = 0.473 fm. 



M^GeV) 


M K 


1 Mi q -2Mj, \ * 
V M K ) 


M K 
M 6 


0.8249 


0.223(3) 


0.296(5) 


0.563(7) 


0.6672 


0.214(6) 


0.267(8) 


0.523(8) 


0.5015 


0.209(8) 


0.252(10) 


0.513(12) 


0.4224 


0.207(12) 


0.239(13) 


0.512(14) 


0.3617 


0.193(15) 


0.247(18) 


0.508(17) 


0.2818 


0.202(22) 


0.234(23) 


0.510(22) 


0.2218 


0.193(26) 


0.232(27) 


0.504(25) 



To verify whether analysis at relatively large quark 
masses would affect the manifestation of the J = 2 state 
and aid to confirm the indication of a resonance, we allow 
the quark mass to be m q > 2m s so that the threshold 
for the decay q 2 q 2 — > {qq){qq) is elevated. The heavy 
quark mass suppresses relativistic effects, which compli- 
cates the interpretation of light-quark states. The re- 
sulting extracted mass ratios are shown in Fig. [5] and 
tabulated in Table UHl 

The behaviour observed for the mass differences be- 
tween the J = 2 and the two-particle states, at large 
quark masses, implies that at larger quark masses, the 
data appear above the two-</> threshold by ~ 95 MeV and 
remains constant in magnitude as the physical regime is 
approached. This trend continues in the physical limit 
where the masses exhibit the opposite behaviour to that 
which would be expected in the presence of binding. 
Again, the positive mass difference could be a signature 
of repulsion in this channel. This suggests that instead of 
a bound state, we appear to be seeing a scattering state 
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FIG. 5: As in Fig. [4] but for larger quark mass. 

in J = 2 channel. Since the mass difference between the 
reported experimental s 2 s 2 mass and the physical 2m^, 
continuum is ~ 20 MeV, the observed signal is too heavy 
to be identified with the empirical /2(1950) or / 2 (2010). 



TABLE III: Hadron mass ratios at larger quark masses. 



M„(GeV) 


U iq -2M 4 
M K 


/Af 4a -2Af^,S* 
^ M K ) 


M K 
A/* 


1.2549 


0.273(2) 


0.346(3) 


0.764(9) 


1.0272 


0.256(4) 


0.308(5) 


0.641(9) 


0.8215 


0.233(6) 


0.277(7) 


0.586(11) 


0.6552 


0.219(7) 


0.269(9) 


0.532(17) 


0.4217 


0.210(9) 


0.243(12) 


0.5182(20) 


0.3625 


0.203(11) 


0.236(14) 


0.5108(24) 


0.2418 


0.194(16) 


0.237(18) 


0.5072(27) 



Using the physical kaon mass, Mk = 503(5) MeV, 
we obtain a mass estimates of 2123(33) (58) MeV and 
2137(39) (64) for the s 2 s 2 tetraquark ground state and 
the second ground-state, respectively. In each case, the 
first error is statistical, and second one is our estimate of 
combined systematic uncertainty including those coming 
from chiral extrapolation and quenching effects. Note 
that we cannot estimate the discretization error since we 
have only one lattice spacing to work with. Given the fact 
that the ratio does not show any scaling violations, we 
could also quote the value of this quantity on our finest 
lattice, which has the smallest error. Nevertheless, order 
2% errors on the finally quoted values are mostly due to 
the chiral extrapolations. The quenching errors might 
be the largest source of uncertainty. Note however, that 



in the case of mass ratios of stable hadrons, this is not 
expected to be very important. It has been shown (3~H 
that with an appropriate definition of scale, the mass 
ratios of stable hadrons are described correctly by the 
quenched approximation on the 1 — 2% level. To this 
end we also calculated the pseudoscalar to vector meson 
ratio Rsp and pseudoscalar to nucleon mass ratio Rsn 
and found that in the physical limit these ratios differ 
about 1% from their corresponding experimental values. 
So we quote our quenching errors to be less than two 
percent. 



IV. SUMMARY AND CONCLUSION 

We presented the results of our investigation on the 
tetraquark systems in improved anisotropic lattice QCD 
in the quenched approximation. The mass of J = 2 state 
was computed using field operators, which arc motivated 
by the non-7T7r and diquark structure. In the quenched 
approximation, our results suggest that our interpola- 
tors have sufficient overlap with fa^ssss) to allow a suc- 
cessful correlation matrix analysis and produced the evi- 
dence that the mass of the lowest-lying state only agrees 
marginally with the mass of /(2010). In the region of 
pion mass which we are able to access, we saw no ev- 
idence of attraction that could be associated with the 
existence of a resonance in J = 2 channel. Since our 
estimated value for the mass of /2(s 2 s 2 ) is marginally 
close to its experimental value, we suspect that might be 
the /(2010) resonance captured by our optimized corre- 
lator. However, on the other hand, our spectral weight 
ratio for two different lattice volumes deviates from one 
(the essential criterion for resonance) with large errors 
for small quark masses, observed state exhibits the ex- 
pected volume dependence in the spectral weight for two 
particles in a box. The ground-state is found to be con- 
sistent with scattering state. Our estimated values serve 
as predictions of lattice QCD in quenched approxima- 
tion. Indeed, our simulation does not include dynamical 
quarks, the final conclusions will have to wait till both 
disconnected correlators and annihilation contributions 
are incorporated. 
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